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Structure and dynamics of the interface between a binary hard-sphere crystal of 

NaCl type and its coexisting binary fluid 
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Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, USA 

(Dated: February 2, 2008) 

Molecular dynamics simulations are performed to study the [100] and [111] orientations of the 
crystal-melt interface between an ordered two-component hard sphere with a NaCl structure and its 
coexisting binary hard-sphere fluid. The diameter ratio of the two types of hard spheres making up 
the mixture is taken to be a = 0.414. This work complements our earlier interface simulations [J. 
Chem. Phys. 116, 3410] for the same diameter ratio at lower pressures where the smaller component 
is immiscible in the solid and the fluid mixture coexists with a pure FCC crystal of large particles. 
Density profiles and diffusion coefficient profiles are presented for the AB interfacial system. We 
find that for this system, the transition from crystal-like to fiuid-like behavior of both the density 
and diffusion constant profiles occurs over a narrower region than that seen in our previous studies 
[J. Chem. Phys. 116, 3410] of the FCC/binary fluid system. But similar to what was found in the 
FCC/binary fluid interface the transition region for the large particle diffusion constant is shifted 
about l.OaA toward the fluid phase relative to that for the small particles. 



I. INTRODUCTION 

The kinetics of crystal growth and nucleation from the 
melt is highly dependent upon the structure, dynarp.- 
ics and thermodynamics of the crystal-melt interfaceo. 
Given the difficulties in obtaining unambiguous infor- 
mation from experiments, most of what is currently 
known about the microscopic phenomenology of such in- 
terfaces is obtained via computer simulationa. Although 
most simulatiMj, aljudies have focused on single compo- 
nent systemsBatl'B'Q'ffl, thece. has. been recent interest in 
multicomponent interfacesE30£3't3. All of these stud- 
ies have involved crystal phases that are either disor- 
dered or pure face-centered cubic (FCC) lattices. In this 
work we present results for the structure and dynamics 
of the interface between an ordered two-component hard- 
sphere crystal with a sodium chloride (NaCl) structure 
and a binary hard-sphere fluid. Such a system can be 
viewed either as as a prototype to understand the inter- 
face between inter-metallic compounds and their coex- 
isting fluid phases or as a model two-component colloidal 
dispersion. 

The hard-sphere interaction was chosen for this study 
since it is an important reference mndel for the study of 
simple liquidslij and liquid mixturcalj. This is especially 
true with regard to phenomena associated with the freez- 
ing transition. For example, it has been recently shown 
that the interfacial free energy of close-packed metals can 
be describedi-ss/ith quantitative accuracy using a hard- 
sphere modeltj. In addition, recent phase boundary cal- 
culations have shown that binary hard spheres form a 
wide range of crystal structures depending on the ratio, 
a = ob/cta, of the diameter of the small spheres (labeled 
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B), (Tb to that of the larger spheres (labeled A), ua- 
A substitutionally disordered FCC crystal is the stable 
phase for 1.0 > a > 0.85 EH while for a < 0.85, only or- 
dered crystal structures are seep-|ta|iie-S|table, including 
AB, ABa and AB13 structuresBBE30E3. A detailed 
study of the disordered FCC crystal /melt interface for 
a — 0.9 has been recently reportedtl. 

In this work, we examine two-component hard-sphere 
mixtures with a considerably larger size asymmetry of 
a — 0.414. This size ratio is significant in the theory of 
binary alloys in that it is the largest asymmetry in which 
the small spheres can be accommodated in the intersti- 
tials of a densest close-packed crystal of larger spheres. 
The phase diagram for this value of a has been deter- 
mined as a function of pressure and mole fraction-using 
MC and MD simulations by Trizac and coworkerstS and 
is reproduced in Fig. |l]. At low pressures, the binary 
fluid coexists with a pure FCC crystal of large spheres, 
whereas at higher pressures (above 50kT/a^) the coex- 
isting solid phase is an ordered 1:1 crystal of the sodium 
chloride (NaCl) type. Earlier cell theory calculations 
also pcedicted the stability of the NaCl at this diameter 



ratioEJ. Other AB structures such as the CsCl and the 
zinc blende have .been shown to be unstable at this di- 
ameter ratioaljcj. Throughout the text we will be using 
AB to also refer to the NaCl structure. 

A detailed study of the low pressure wre FCC/binary 
fluid system has been recently reportedllj. In that work, 
the coexistence between the crystal and an approxi- 
mately 1:1 binary mixture was examined, correspond- 
ing to a pressure of 20.1fcr/(T^, which is approximately 
twice the coexistence pressure {11.55kT/a) of the pure 
single component system. As the pressure is increased 
the mole fraction of large spheres in the fluid phase, X-^, 
decreases, and at a pressure of about bOkT/a^ the fluid 
coexists with an AB crystal of NaCl type (see Fig. |l|). In 
this work, we examine in detail the microscopic structure 
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FIG. 1: Pressure-concentration phase diagram of the bi- 
narv-hard-sphere system with a = 0.414. [Reprinted from 
Ref.LJ by permission of the pubhsher, Taylor and Francis, Ltd 
(www.tandf.co.uk/journals)]. Note that to make the phase 
coexistence lines easier to distinguish, the pressure is plotted 
against the cube root of the large sphere mole fraction. The 
phase point of the present study is shown as a filled cirqki. 
The open circle shows the conditions of our previous studyt3. 



and dynamics of the interface between the high pressure 
AB crystal and its melt. To do this, we have chosen 
a point in the phase diagram with xj^ ~ 0.097 . At 
this mole fraction, the fluid coexists with the NaCl crys- 
tal at a pressure of 53kT/a\. Details of the simulation 
methodology and interface equilibration procedure are 
outlined in the next section followed by presentation of 
the results of the study in Section III and concluding 
remarks in section IV. 



II. SIMULATION DETAILS 

Molecular-dynamics simulations are performed on a 
two-component system of hard spheres of differing di- 
ameters, a A and as, with a a < <Jb (Type A particles 
are assumed to be the larger spheres). The interaction 
potential between two spheres is defined by the following 
pairwise interaction potential 



oo 




r < a. I 
r > a I 



(1) 



where i,j e {A,B}, r is the distance between the cen- 
ters of the two interacting spheres, and (Jij = {(7i + (jj)l2 
is the distance of closest possible approach between two 
spheres with diameter ai and tTj, respectively. The sys- 
tem is completely defined by specifying the diameter ra- 
tio a = (Ts/ca, the mole fraction of large particles, xa 



and the total number density, p. The total volume occu- 
pied by the hard spheres relative to the volume available 
to the system is given by the packing fraction. 
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-p[xA + (1 - XA)a^ 



(2) 



where p = pA + Pb — N/V is the total density, xa is 
the mole fraction of the larger species, and a — as/o'A 
is the diameter ratio. 

The procedures for interface construction and equili- 
bration of binary interfaces employed in Xhis study are 
similar to those used in our earlier workli3 on the low- 
pressure coexistence in this system (single component 
FCC {xa = l)/binary fluid mixture). Since the gen- 
eral interface preparation process is described at length 
in that work, only those details specific to the current 
system are described here. Tlie interested reader is en- 
couraged to consult referencai3 for a more complete de- 
scription. 

To construct an interface, fluid and crystal blocks are 
prepared separately using the calculated coexistence con- 
ditions. At a pressure of 53 kT / a\, we independently de- 
termined the packing fractions of the coexisting crystal 
and fluid to be t]'^ = 0.683 and rj^ — 0.490, respectively 
In the preparation of a fluid block, it is a usual prac- 
tice to initially position particles in a lattice at a den- 
sity lower than coexistence. As the system is allowed to 
equilibrate the lattice melts, giving a fluid conflguration, 
which is then compressed to the coexistence density. For 
the system under study here this procedure is not feasi- 
ble since the size asymmetry {a ~ 0.414) and the mole 
fraction of small particles (A^ = 0.903) in the fluid sys- 
tem are both large. (It should be noted that for this 
value of a the large sphere volume is over 14 times that 
of the small particle.) To construct the fluid phase we 
began with an FCC lattice of small particles at a number 
density equal to the desired total number density of the 
fluid mixture. A number of particles, corresponding to 
the target mole fraction of large spheres, are then chosen 
at random from this lattice. A molecular dynamics run 
is started from this initial conflguration and the diame- 
ter of randomly chosen particles is periodically increased 
until the correct large sphere diameter is attained. The 
amount of increase in the diameter at each stage depends 
on the maximum increase that is possible without cre- 
ating particle overlap. The preparation of the crystal is 
straightforward as the small particles easily inserted into 
the interstitial sites of the large sphere FCC lattice. 

After equilibration of separate crystal and fluid sys- 
tems, they are placed in contact within the simulation 
box. Due to the periodic boundary conditions, two in- 
terfaces are formed. Note that packing fraction used in 
the preparation of the initial fluid block is slightly dif- 
ferent from the predicted rjf because a gap of Icr^ is 
placed between the crystal and fluid blocks to avoid any 
initial overlap that may occur when the two blocks are 
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TABLE I: Number of particles and dimensions of the 
simulation box 







AT/ 


Lx/UA 


Lyja A 


Lz/(TA 


[100] 


7056 


14976 


10.41 


10.41 


53.02 


[111] 


7200 


15552 


10.51 


10.92 


51.44 



combined. The fluid particles are then allowed to move 
while the large particles are fixed so the initial gaps are 
then filled with the fluid particles. The initial fluid pack- 
ing fraction is adjusted until an unstressed bulk crystal 
is obtained when the two blocks are combined and equi- 
librated 

In this study, we examine both the [100] and [111] 
orientations of this NaCl crystal/binary fluid interface. 
For reference, we define the 2;-axis to be perpendicular 
to the interfacial plane. Periodic boundary conditions 
are applied in x, and z directions. The length in the 
^-direction, L^, is a sum of lengths of the two separate 
blocks in the z-direction and the 2cryi-gap initially left 
between the two blocks. The total number of particles 
used are 22,032 and 22,752 for the [100] and [111] ori- 
entations, respectively. The details of the system sizes 
used are summarized in Table As the systems stud- 
ied ar|e,large, we have implemented the cell method tech- 
niqueEj to efficiently carry out the molecular-dynamics 
simulationscj. Equilibration was done for about 4 x 10^ 
collisions per particle (cpp). During sampling, the run 
was divided into blocks of 2600 cpp each. The sampling 
run was over a length of 30 blocks, but since each block 
gives two independent measures of interfacial properties 
when the system is properly folded at the center of the 
crystal, the results reported are averages taken over 60 
samples. 

In the course of such simulations the position of the 
interfacial plane can shift due to Brownian motion or 
due to melting or freezing of the crystal. To prevent 
such motion from artificially broadening the interfacial 
profiles it is necessary to monitor the position of the in- 
terfacial plane during the sampling runs. The standard 
measure of interfacial position for a planar interface is 
the Gibbs dividing surfaceEl defined for a multicompo- 
nent system as the value of z where the surfaces excess 
number of "solvent" particles is zero. However, accurate 
calculation of the Gibbs dividing surface requires rela- 
tively long simulation runs and is then unsuitable for a 
method of monitoring the time dependence of interfacial 
position. In our previous studieg3't3, we find that a more 
suitable measure can be obtained from measurement of 
orientational order parameter profile as a function of z. 
The orientational order parameter is defined as follows 



where n = 4 for the [100] orientation and n — Q for 
the [111], j, j and k are nearest neighbor large particles, 
Oxy(i,j,k) is the bond angle formed by i,j and k pro- 
jected on the x,y plane, and Nz is the total number of 
particles that form bond angles. The average is taken 
over the number of angles found between z — Az/2 and 
z + Az/2, where Az is equal to the layer spacing of the 
bulk crystal. 

We show in Fig. ^ the order parameters, and gg, 
of the two particle types for the [100] (upper panel) and 
[111] (lower panel) orientations. As expected 94 (gg) is 
small in the [111] ([100]) interface where 6-fold (4-fold) 
symmetry dominates. We define the interfacial position 
relative to the midpoint of the orientational order pro- 
file for the large particles. That is, z = in all of the 
z— dependent profiles presented in this study is defined 
as the point at which the orientational order parameter 
has decayed halfway from its crystal to fluid value. This 
order parameter is suitable as a measure of interfacial po- 
sition since it is smoothly monotonically decreasing and 
can be calculated accurately for very short runs. The pa- 
rameter profiles of the small particles are not smoothly 
varying because at the interfacial region, some number 
of small particles cluster together to occupy large parti- 
cle vacancies at the interfacial region (as will be seen in 
the density plots presented in the next section) disrupt- 
ing the smooth transition from crystal-like to fluid-like 
value of the orientational order parameter. 

Analysis of the interfacial position as a function of 
time shows that during the equilibration run the crystal 
exhibits some initial growth, but quickly stabilizes before 
the averaging runs are begun. Brownian motion of the 
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FIG. 2: Large (dotted) and small (solid) particle orientational 
order parameter profiles, 54 (circle) and qa (square) for the 
[100] and [111] interfacial orientations. The point z = is 
the location of the interfacial plane as calculated from the 
large particle order parameter profile. 
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FIG. 3: Large (dotted) and small (solid) particle density pro- FIG. 4: Filtered density profiles for the large (dotted), small 
files of the NaGl/binary fluid interface for the [100] and [111] (solid) and total (dashed) densities in the (a) [100] and (b) 
orientations. The distance between vertical dotted lines is [HI] interfacial orientations, 
equal to the crystal layer spacing in [100] and twice the spac- 
ing in [111] 



solid phase, as monitored by motion of the inner layers 
of the crystal, was found to be neghgible due to the large 
system size and no correction was necessary. 



III. SIMULATION RESULTS FOR THE [100] 
AND [111] INTERFACES 

A. Structure: Density profiles and contour plots 

The structural variation of the system across the in- 
terface is determined by calculating the density profile 
for each particle type. 



<N,iz)> 

LxLyl^Z 



(4) 



where i denotes a particle type, Az is 1/25 of the crystal 
layer spacing, <Ni{z)> is the average number of particles 
of type i in the region between z ~ Az/2 and z + Az/2. 

The density profiles of the two particle types are shown 
in Fig. ^(a) and (b) for the [100] and [111] orientations, 
respectively. As expected for an NaCl-type lattice, the 
small and large particle peaks are in phase in [100] and 
exactly out of phase in [111]. This is similar to the reg- 
istry of particle density peaks found at the interfacial 
regie* of the lower pressure pure FCC/binary fluid sys- 
temt3 . Due to the higher pressures in this study, the 
crystal peaks are much sharper than those seen in the 
lower pressure, binary systemllHI or in the single compo- 
nent interfaceQ. The small side peaks in the NaCl den- 
sity are due to the filling of large particle vacancies in the 



lattice structure with several smaller particles - as dis- 
cussed below. In order to reveal any change in the lattice 
spacing through the interface and to index the interfacial 
planes for later use, vertical dotted lines separated by the 
bulk crystal lattice spacing were added to Fig. ||. Labels 
a to i in the [100] profile marks some layers whose cross- 
sectional density distributions have been determined and 
will be discussed later. For both orientations, there is no 
discernible change in the lattice spacing as the interface 
is traversed from crystal to fluid, in contrast to what was 
seen in our previous studies where Jie-crystal was either 
a pure or disordered FCC latticeCl'LilEj. In those studies 
there was a significant increase in the lattice spacing in 
the [100] orientation as the fiuid side of the interface was 
approached. 

The oscillations in the fine scale density profiles shown 
in Fig. ^ make it difficult to see the overall trend in bulk 
density, so we have processed these profiles using a Finite 
Impulse Response filtcitlc3 to reveal the non-oscillatory 
component of the density variation. The resulting fil- 
tered density profiles are shown in Fig. ^. The 10-90 
width of these bulk density profiles provides a measure 
of the interfacial width. (The 10-90 width of a mono- 
tonically varying interfacial profile is the distance over 
which the profile changes from 10% to 90% of the higher 
of two coexisting bulk values, relative to the lower bulk 
value.) The 10-90 width derived from the large particle 
density profile of the [100] orientation [see Fig. ^(a)] is 
2.6(Tyi, corresponding to the region between z = —2.2a a 
to z — 0.4(7^. For the small particles the 10-90 width 
is smaller at 2.3crA and the 10-90 region [z = — I.Tcta 
to z = O.Gcta) is shifted slightly toward the fluid, rela- 
tive to the 10-90 region of the large spheres. Combining 
these two regions, the interfacial region of the [100] ori- 
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FIG. 5: Large particle density contour plots parallel to the 
interfacial plane for different layers of the [100] interface. The 
layers are as labeled in Fig. H(a). 



FIG. 6: Small particle density contour plots parallel to the 
interfacial plane for different layers of the [100] interface. The 
layers are as labeled in Fig. ^(a). 



entation defined by the densities has a width of 2.8a-A- 
The total interfacial width defined for the densities for 
the [111] orientation (2.9cr^) is not significantly differ- 
ent than that for [100]. The interfacial widths of the 
large particle density profiles are identical to those found 
in ou&-carlicr lower pressure FCC/binary fiuid interface 
studyll3 for the same diameter ratio, but narrower than 
the 3.3(7 widths found in the single component systemQ. 
In contrast, the small particle interfacial widths found 
here are much smaller than those found in the lower pres- 
sure study, where they were found to be about 3.3a a- As 
a consequence, the overall interface for the NaCl/fluid is 
slightly narrower than the lower pressure pure FCC/fiuid 
interface. 

To get a detailed understanding of the structural 
transition across the interface between the two coexist- 
ing phases, we have determined the cross-sectional den- 
sity distributions within layers parallel to the interfacial 
plane. For each layer, the 2-D density distribution is 
defined as 



<Nf{x,y)> 
AxAyAz 



(5) 



where i denotes a particle type. Ax = Ay = 0.12ct^, Az 
is the crystal layer spacing, which is 0.74a^ for [100]. 
The average number of particles of type i in the volume 
given by AxAyAz is < N--{x, y) >■ In Figs. H] and ^ are 
density contour plots of the [100] interface for the large 
and small particles, respectively. Layers a to i are as la- 
beled in Fig. S(a), where a is deep into the bulk crystal. 



i is in the bulk fluid and b to h are interfacial regions. 
The decrease in density peak height of the large (Type 
A) spheres in Fig. ^ starting at layer b is initially due 
to the presence of lattice vacancies as shown in Fig. ^. 
Small particle vacancies also start to appear in layer b of 
Fig. H. We also find by comparing Figs. ^ and || for layers 
c and d that the side peaks appearing at these layers are 
due to the accumulation of small particles in the large 
sphere vacancies. The structure of this vacancy filling is 
interesting in that the large particle is typically replaced 




FIG. 7: Snapshot of a portion of layer c as labeled in Fig. ^(a) 
showing small particles occupying a large particle vacancy. 
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by 6 small particles (although a small number of vacan- 
cies filled with 5 or 7 small spheres do occur) with little 
disturbance to the surrounding lattice. This can be seen 
in Fig. where a snapshot of one of these filled vacancies 
in layer c is shown. A uniform 2-D density distribution 
begins to develop in layer G for both particle types indi- 
cating that, although the z— dependent density profiles 
still has some oscillations in this region, the structure is 
that of an inhomogeneous fluid at a wall. 

Of particular interest to materials scientists is the de- 
gree of interfacial segregation - the preferential adsorp- 
tion (or desorption) of one component (usually the "so- 
lute" ) at the interface. This quantity is defined relative 
to the Gibbs dividing surface. TheJGIibbs dividing sur- 
face of a planar interface is definedEl as the plane along 
the z— axis giving a vanishing surface excess solvent par- 
ticle number, Tsolvent' defined in the equation 

^solvent = l^ent^ + l^^nt(i. -z) + V^^^^^^^ 

(6) 

where TV^'^^^®'^^ the total number of solvent particles 
spheres, A is the area of the interface, ^solvent 
^solvent g^j.g |.]^g bulk densities, z is the location of the 
interface assuming the length of the simulation box runs 
from to Lz- Defining the small particles as the 'sol- 
vent', we find that the Gibbs dividing surfaces are at 
z = —0.49(7^ and z = — 0.93(t^ relative to the position 
calculated from the large particle order parameter for 
the [too] and [fff] orientations, respectively. Surface 
excess of the 'solute' particles, was found to be neg- 
ligible, indicating the absence of interfacial segregation, 
a result that is consistent fQtxi|thcr crystal/melt systems 
that have been investigatedEruO. 

B. Transport: Diffusion coefficient profiles 

Inhomogeneities in the transport properties within 
the interfacial region can be examined by calculating 
z— dependent diffusion coefficient profiles, defined for a 
particle of type i by 



0.06 



D,{z) = lim 



f 



t^oo m^iz) dt 



^ ([r,(i)-r,(io)P). (7) 



The term in the summation is the mean-squared dis- 
placement over a time interval t — to of a total of Ni type 
i particles located between z — Az/2 and z + Az/2 at 
time toi where Az is the layer spacing in [100] and is 
twice the layer spacing in [Iff]. 

Figure ||(a) shows the diffusion coefficient profiles for 
the current study (dotted lines), including for compari- 
son the results previously r-eeorted for the lower pressure 
FCC/binary fluid interfacey (solid fines). Only the [100] 
results are shown as the diffusion profiles for the [Iff] in- 
terfaces are not statistically different. The error bars are 




FIG. 8: (a) Diffusion coefficient profiles for the [100] orien- 
tation of the FCC/ffuid (dashed) and the AB/fluid (solid) 
interfaces for both particle types (circle for large and square 
for small); (b) Corresponding diffusion coefficient plots scaled 
to be zero in bulk crystal and unity in bulk fluid. 



small and so are not shown for clarity of the plots. The 
bulk fluid value for the large particles in the lower pres- 
sure FCC/binary fluid system is 

0.012(A:TcryTO)i/2 and 
that for the small particles is 0.050(A;T(t^/to)^/^. Since 
the AB/fiuid system has a higher pressure and larger 
fluid packing fraction, the bulk fluid diffusion coefficient 
values are lower: 0.008(A:Tcr^/m)^/'^ for the large parti- 
cles and 0.029{kTa\/my^^ for the smaU particles. The 
difference in magnitude between the small and large par- 
ticle diffusion constants makes it difficult to compare 
the two diffusion profiles. For a clearer comparison in 
Fig. ^(b) we plot the data in Fig.^(a) normalized rela- 
tive to the bulk fluid values. Traversing the system from 
fluid to crystal, we find a region of width greater than 
laA where the small particles have nonzero diffusion co- 
efficient while the large particles have have effectively 
zero diffusion. Both high- and low-pressure systems ex- 
hibit this shift in the change from crystal-like to fluid-like 
motion of the two particle types. 

As was done for structural transition, we can also de- 
fine the extent of dynamical transition by determining 
the 10-90 region from diffusion coefficient profiles. From 
the diffusion coefficient profile of the large particles this 
region starts from z — — 0.7cr^ up to z = l.Scr^, result- 
ing to a width equal to 2.0aA and centered at z = O.SaA- 
The small particles define an interfacial region that starts 
from z = —2.2aA and ends at z = O.QaA- These bound- 
aries give a width of 2. 8a a, which is 40% greater than 
the width from diffusion of the large particles. Also, the 
midpoint is shifted by l.lcr^ to the fiuid side relative to 
the midpoint of the region defined by the large particles. 

As in our previous studies on hard-sphere 
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FIG. 9: Filtered density, diflFusion and orientational order 
parameter profiles for the [100] interface. All profiles are 
scaled such that they go from unity in the crystal to zero in 
the fluid phase. 



interfacesilii the regions defined each of the indi- 
vidual density and diffusion coefficient profiles do not 
all coincide, so we take a union of these 10-90 regions 
to define the full interfacial region. Taking into account 
all profiles, the NaCl/binary fluid hard-sphere interface 
studied here has an overall 10-90 width of 3.5(Ta, 
corresponding to layers c to g in Fig. |^. As perhaps 
expected for the higher pressure system, this interfacial 
region is narrower than the pac found for lower pressure 
FCC/binary fluid interfaced, which was found to be 
A.SaA- To summarize all of the profile data for quick 
inspection we show in Fig. ^ all order parameters 
profiles, normalized such that we have all values equal 
to unity in the bulk crystal and zero in the bulk fluid for 
the [100] interface. Except for differences in the overall 
interfacial width, this plot is qualitatively quite similar 
to the corresponding plot for our earlier low-pressure 
FCC/binary fluid simulations in that the transition of 
densities for both particle types and the diffusion of 
the small particles is observed over approximately the 
same region, while the transition for the large particle 
diffusion is shifted by about laA (l-3crA in the earlier 
study). Another notable feature in this plot is that the 



transition for the orientational order parameter, which 
we use to locate the interfacial plane, occurs at about 
the center of these two transition regions. 



IV. SUMMARY 

Using molecular-dynamics simulation, we have inves- 
tigated the structure and dynamics of the [100] and [111] 
crystal-melt interfaces of the AB(NaCl) /binary fluid sys- 
tem for a two-component hard-sphere system in which 
the ratio of small to large particle diameter is a = 0.414. 
This system was at a pressure of 53kT/a\ that is at the 
lower range of the NaCl crystal/fluid coexistence-region. 
These simulations complement our earlier worUl3 on the 
pure FCC/binary fluid interface found in this same sys- 
tem at lower pressure {20.1kT/a\). We find that the 
higher pressure AB /binary-fiuid interface has a narrower 
interfacial region of 3.5aA compared to the lower pres- 
sure FCC/binary fluid system at the same diameter ra- 
tio, which had an interfacial region of width 4.8(t^. In 
addition, the crystal side of the higher pressure binary in- 
terface exhibited much higher vacancy defect concentra- 
tions than either the low pressure binary system or the 
single-component interfaceQ. In the interfacial region, all 
vacancies in the large particle lattice were found, with lit- 
tle distortion in the surrounding lattice, to be flUed with 
an average of 6 small particles. Similar to what was 
seen in the lower pressure FCC/binary fluid interface, 
the transition regions for both density proflles and the 
small particle diffusion constant are approximately coin- 
cident whereas the diffusion profile for the large particles 
is shifted relative to the others by about Icta toward the 
liquid side of the interface. 
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